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ABSTRACT 


A  theoretical  investigation  (§£  initio)  and  experimental  data  have  been 
combined  in  an  attempt  to  characterize  the  rhodium  carbonyl  hydride  species 
which  has  been  shown  previously  in  these  laboratories  and  in  those  of  Solymosi 
to  be  a  long-lived  surface  species  for  supported  rhodium  in  the  catalytic 
hydrogenation  reactions  for  00  and  (X^.  The  calculations  have  predicted  that 
a  bent  rhodium  carbonyl  hydride  species  is  much  more  stable  than  a  rhodium 
formyl  species.  They  have  also  predicted  that  the  Rh-H  stretching  frequency 
should  be  more  than  300  an"’1  lower  than  the  C-0  stretching  frequency  for  the 
rhodium  carbonyl  hydride  species.  This  prediction  is  at  odds  with  assigned 
Rh-H  stretching  frequencies  for  Rh  organometallic  complexes  which  place  Rh-H 
and  C-0  stretching  frequencies  in  the  same  region  of  the  spectrum  (2000  -  2200 
cm-1) .  The  experimental  data  are  not  conclusive  on  this  point.  The 
calculations  do  satisfactorily  rationalize  the  enhanced  dissociation  of  the 
C-0  bond  for  the  rhodium  carbonyl  hydride  species  caused  by  the  presence  of 
the  hydride  moiety. 
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Work  in  these  laboratories  in  recent  years  has  focused  on  the  search  for 
intermediates  in  the  reactions  of  hydrogen  with  carbon  monoxide  '  and  carbon 
dioxide5-5  over  supported  rhodium  films;  infrared  spectroscopy  has  been  the 
primary  analytical  probe  enployed  in  these  studies.  A  key  Rh  surface  species 


which  has  been  identified  by  infrared  for  both  of  the  hydrogenation  reactions 
is  the  carbonyl  hydride  shown  below  which  exhibits  an  infrared  band  in  the 
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Rh 


region  2020-2045  cm-1  dependent  upon  the  nature  of  the  support  and  whether  one 

or  two  hydride  surface  atoms  are  present.5  Solymosi  and  coworkers  were  the 

first  to  postulate  the  existence  of  the  rhodium  carbonyl  hydride  species 

6  7 

produced  during  the  reduction  of  CX>2  over  supported  Rh.  '  It  has  been 
concluded  that  the  presence  of  the  hydride  moiety  causes  a  shift  of  the 
frequency  for  the  00  stretching  mode  to  a  lower  value  (2020-2045  cm  1)  than 
that  for  the  "linear"  Rh  -CO  species  (2040-2080  cm-1)  due  to  electron  donation 

A 
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into  the  ff*  orbital  of  the  CO  moiety.  This  satisfactorily  explains  the 
resulting  weakened  CO  bond,  leading  to  its  facile  dissociation  when  hydrogen 

r 

is  present.  This  argument  was  challenged  by  Iizuka  and  coworkers  who 

concluded  that  the  low  frequency  CO  mode  was  not  due  to  a  carbonyl  hydride 

species,  but  rather  merely  to  the  usual  linear  Rh-CO  species  at  low 

coverage.1-®' 11  However,  the  argument  of  Iizuka  and  coworkers  seemed 

improbable  to  us  because  the  linear  Rh-CO  species  is  always  accompanied  by  a 

bridged  Ri^CO  species  having  a  broad  infrared  band  in  the  1840-1920  cm  1 

12 

region;  both  species  are  formed  on  metallic  cluster  sites.  The  only  surface 
species  detected  by  infrared  for  the  hydrogenation  of  (X>2  over  supported  Rh  is 
the  one  (carbonyl  hydride)  exhibiting  an  infrared  band  in  the  2020-2045  cm  1 


region.  This  was  borne  out  by  deuterium  experiments  in  these  laboratories 
which  showed  conclusively  (a  small  observed  red  shift)  that  the  species 

producing  the  infrared  band  in  the  2020-2045  cm  ^  region  was  indeed  a  rhodium 

4  / 

carbonyl  hydride.  In  a  very  recent  study  Solymosi  and  Pasztor  have  shown 

that  the  presence  of  the  rhodium  carbonyl  hydride  species  during  the 

hydrogenation  of  CX>2  prevents  the  disruption  of  Rh  clusters  such  that  the 

usual  gem-dicarbonyl  species  always  detected  during  CO  hydrogenation  on 

+  13 

isolated  Rh  sites  can  not  form. 

It  should  be  noted  that  a  formyl  species  might  also  be  considered 


Rh 


as  a  possible  intermediate  in  the  hydrogenation  of  CO  and/or  C02.  However, 

this  surface  species  has  not  been  detected  for  supported  Rh  catalysts . 

Indeed,  Yates  and  Cavanagh  have  reported  a  Herculean  effort  in  which  the 

14 

formyl  species  was  sought,  but  could  not  be  observed.  One  might  pose  the 
question;  why  is  the  rhodium  carbonyl  hydride  species  observed,  while  a 
rhodium  formyl  species  is  not?  The  purpose  of  the  current  work  was  to  address 
this  question  using  high  quality  afe  initio  computations  and  Fourier  transform 

r 

infrared  spectroscopy .  This  work  represents  one  of  the  initial  attempts  at 
using  ab  initio  molecular  orbital  calculations  for  transition  metal  species 
relevant  to  surface  chemistry. 

THEORETICAL  AND  EXPERIMENTAL  SECTION 

For  the  computations  in  this  study,  the  36  core  electrons  of  Rh  were 
replaced  by  a  relativistic  effective  core  potential  (RECP)  as  generated  by  Hay 
and  Wadt.^  The  (5s4d5p)  valence  orbitals  were  represented  by  a  (3s3p4d) 
primitive  Gaussian  basis  contracted  to  (lslpld)  to  form  a  basis  set  denoted 


BAS1.  A  basis  set  denoted  BAS2  was  formed  by  splitting  off  the  outermost 

17 

primitive  to  form  a  (2s2p2d)  contraction.  A  standard  6-31G  basis  was  used 
for  C  and  0  in  BAS1,  while  a  set  of  diffuse  s  and  p  functions  was  added  to  C 
and  0  (6-31+G)  in  BAS2.  Hydrogen  was  described  by  four  primitive  gaussians; 
the  inner  three  were  contracted  as  found  in  the  6-31G  basis,  while  the  outer 
gaussian  was  varied  independently  with  an  exponent  of  0.1219.  The  basis  set 
BAS1  was  used  to  calculate  optimized  geometries  and  force  constants.  All 
geometric  variables  were  optimized  within  the  assumed  symmetry  point  group, 
and  frequencies  were  computed  by  taking  finite  differences  of  analytical  first 
derivatives.  The  basis  set  BAS2  was  used  to  determine  electron  populations. 
The  "best”  energies  to  be  discussed  were  determined  from  MP3/BAS2  energies  and 
zero  point  correction  from  vibrational  frequencies  (corrected  by  a  factor  of 
0.89) 19  at  the  BAS1  level. 

The  experimental  procedures  employed  in  these  laboratories  for  infrared 
investigations  of  supported  Rh  catalyst  films  have  been  discussed  in 
detail.1  ^  Briefly,  sanples  consisting  of  thin  films  (4.3  mg  cm  )  of  2.2% 
Rh/TiC>2  were  prepared  by  spraying  a  slurry  of  RhCL^Sl^O  and  Ti02  (Degussa 
Titanium  Dioxide  P25,  50  m  g  )  in  water/acetone  from  a  specially  designed 

r 

atomizer  onto  a  25  mm  CaF2  infrared  window  held  at  ca.  360  K.  The  solvents 

rapidly  vaporized  leaving  a  film  of  RhClg* 

window.  The  window  containing  the  catalyst  was  mounted  in  a  specially 
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designed  Pyrex  infrared  cell  which  was  evacuated  to  ca.  10  torr  using  an 
ion  pump.  The  sample  was  then  reduced  with  85  torr  doses  of  hydrogen  gas  at 
480  K  for  periods  of  5,  5,  10,  and  20  min,  each  reduction  cycle  followed  by 
evacuation.  After  heating  for  1  h  at  490  K  under  vacuum,  the  sample  was 
evacuated  overnight  to  ca.  10  torr  at  ambient  tenperature.  Then  a  77.4  torr 
mixture  of  C02  and  H2  at  a  ratio  of  1:4  was  added  to  the  cell,  and  the 


3H20/Ti02  adhered  tightly  to  the 


methanation  reaction  was  initiated  at  440  K  for  1  h  followed  by  evacuation  to 

_3 

10  torr  at  ambient  temperature.  At  this  point  spectrum  a  in  Figure  2  was 
recorded.  Then  100  torr  D2  was  dosed  into  the  cell,  left  overnight,  at 
ambient  temperature,  and  spectrum  b  in  Figure  2  was  recorded.  These  spectra 
were  recorded  on  an  IBM  IR/32  Fourier  transform  infrared  spectrometer  (500 
scans)  at  4  an-1  resolution  using  a  DTGS  detector. 

RESULTS 

The  computed  optimized  geometries  of  the  various  Rh  species  and  CO  and 

HCO  are  given  in  Figure  1.  Molecular  and  electronic  symmetries,  total 

energies,  and  zero  point  energies  for  the  relevant  species  are  presented  in 

Table  I.  Table  II  summarizes  the  relative  energies  of  the  various  species  on 

the  HRhCO  potential  energy  surface.  It  is  clear  that  the  bent  form  (86.6°)  of 

the  rhodium  carbonyl  hydride  was  predicted  to  be  considerably  more  stable  than 

2 

the  linear  form  or  the  rhodium  formyl  species.  The  D  state  for  Rh  could  not 
be  determined  at  the  BAS2  level  because  of  the  collapse  of  the  wave function  to 
one  with  extensive  spin  contamination.  Table  III  presents  computed  bond 
energies  for  the  species  relative  to  fragments  in  their  most  stable  electronic 
states.  Table  IV  provides  calculated  vibrational  frequencies  and  intensities 
for  the  Rh  species  and  CO  and  HCO;  the  vibrational  frequencies  for  RhCO  could 
not  be  determined  due  to  convergence  problems. 

Figure  2  shows  the  infrared  spectra  for  the  rhodium  carbonyl  hydride  and 
rhodium  carbonyl  deuteride  species  adsorbed  on  a  2.2%  Rh/Ti02  catalytic  film. 

DISCUSSION 

The  ab  initio  computations  for  the  various  rhodium  species  reported  in 
this  work,  of  course,  refer  to  gas-phase  species,  rather  than  true  surface 


species.  Thus  variation  of  the  calculated  energies,  bond  energies,  and 
vibrational  frequencies  from  the  true  values  for  the  adsorbed  species  are  to 
be  expected.  Nevertheless,  it  is  believed  that  the  theoretical  data  generated 
in  this  work  can  provide  at  least  a  qualitative  description  of  the  potential 
energy  surface  for  the  rhodium  carbonyl  hydride  and  rhodium  formyl  species. 

Some  experimental  data  do  exist  which  can  be  compared  with  the  ab  initio 

predictions.  Using  mass  spectral  data,  and  based  upon  an  empirical 

relationship  between  electron  affinities  and  bond  dissociation  energies  of 

diatomic  metal  hydrides,  Squires  has  estimated  the  Rh-H  bond  energy  to  be  53.7 
22 

kcal/mol.  At  the  MP3/BAS2  +  ZPC/BAS1  level  of  theory,  which  was  the  best 

level  employed  in  this  work,  the  Rh-H  bond  energy  is  predicted  to  be  73.1 

kcal/mol  (Table  III) .  However,  the  computed  *Rh-H  bond  energy  which  was  48.5 

kcal/mol  is  in  better  accord  with  the  experimental  estimate.  Even  though 

3Rh-H  is  the  ground  state  for  gaseous  Rh-H,  the  1Rh-H  may  be  a  more  realistic 

23 

model  for  the  supported  catalyst.  Connor  ££  al.  concluded  that  the  strength 
of  the  Rh-CO  bond  in  Rhg (CO) ^  was  39.6  kcal/mol,  and  Somorjai  et  al.^ 
estimated  the  energy  of  the  Rh-CO  bond  for  Rh(100)  and  Rh(lll)  single  crystals 
(from  thermal  desorption  studies)  to  be  31  and  29  kcal/mol,  respectively. 

Table  III  shows  that  the  Rh-CO  bond  energy  was  predicted  to  be  11.0  kcal/mol 

4 

for  Rh-CO  and  28.6  kcal/mol  for  the  carbonyl  hydride  species  at  the  best 

level  of  theory  employed  here.  A  recent  X-ray  structure  determination  has 

provided  a  value  of  1.56  A  for  the  Rh-H  bond  distances  in  the  complex 

RhfH^ItNHfSiMSgCHgd-Pr^^]  w^ch  c^-ose  to  the  value  (1.547  A)  predicted 

by  the  sb  initio  calculations  for  the  bent  carbonyl  hydride  species.  Typical 
25 

Rh-C  bond  lengths  fall  in  the  range  of  1.78  to  2.01  A  which  indicates  that 
the  jfc  initio  procedure  does  well  at  predicting  Rh-CO  bond  lengths  (see  Figure 
1) .  The  absorption  frequencies  calculated  (Table  IV)  for  CO  and  HCO  were 


reasonable  for  an  SCF  M3  method;  the  bending  mode  was  particularly  well 
reproduced  for  both  HCO  and  DCO. 

It  is  interesting  that  the  jab  initio  calculations  predict  the  Rh-H  bond 
to  be  very  strong  (92.9  kcal/mol)  and  the  Rh-H  stretching  mode  to  exhibit  the 
most  intense  band  in  the  infrared  spectrum  of  the  carbonyl  hydride  species. 

In  fact,  the  presence  of  CO  is  predicted  to  enhance  the  strength  of  the  Rh-H 
bond.  Repeated  attenpts  have  been  made  in  these  laboratories  and  elsewhere  to 
observe  the  Rh-H  stretching  frequency  for  supported  Rh  catalysts.  If  a 
Rh/A^Og  or  Rh/Ti02  film  is  prepared  in  the  manner  discussed  in  the 
Theoretical  and  Experimental  Section,  and  then  exposed  to  ^  or  D2  at 
pressures  up  to  100  Torr,  with  no  CO  or  C02  present,  no  band  which  can  be 
obviously  attributed  to  Rh-H  or  Rh-D  is  observed.  In  fact,  there  is  no  band 
at  all  in  the  1800-2200  cm  1  region.  Yet  this  does  not  preclude  the  presence 
of  CO  as  a  stabilizing  influence  on  a  Rh  hydride,  and  many  workers  have 
assigned  bands  in  the  2000-2200  cm  region  for  organometallic  rhodium  hydride 
complexes  to  the  Rh-H  stretching  mode.  Also,  Iwasawa  and  Sato^  claim  to 
have  detected  a  "mononuclear,  rhodium-monohydride  species"  with  an  associated 
infrared  band  at  2048  cm-1  during  the  hydrogenation  of  Rh(r^-C3H5)3/Ti02. 

r 

The  frequencies  predicted  for  the  Rh-H  and  Rh-D  stretching  modes  of  the 
carbonyl  hydride  (deuteride)  species  led  us  to  reexamine  the  infrared  spectra 
of  these  species;  previously  only  one  band  in  the  2020-2045  cm-1  region  had 
been  noted  and  assigned  to  the  C-0  stretch  for  the  species.1-9  The  data  in 
Figure  2  clearly  show  that  the  band  at  2037  cm-1  which  shifts  to  2032  on-1 
upon  deuterium  exchange  must  be  assigned  to  a  CO  stretching  mode,  because  if 
it  were  due  to  an  Rh-H  stretch,  the  isotopic  shift  would  be  ca.  600  cm  1  to 
lower  frequency.  However,  a  broad  band  is  observed  in  spectrum  b  of  Figure  2 


—1  “1 
at  ca.  1450  an  which  is  more  intense  than  the  wing  of  the  intense  1615  cm 

band  in  spectrum  a.  This  could  possibly  be  due  to  a  shift  of  a  weak  component 

of  the  2037  an-1  band  noted  in  spectrum  a  of  Figure  2.  In  other  words,  the 

bands  for  the  C-0  and  Rh-H  stretching  modes  of  the  carbonyl  hydride  species 

may  fortuitously  overlap.  However,  the  ab  Initio  calculations  predict  that 

the  Rh-H  stretching  frequency  should  be  near  1643  an  ,  while  the  C-0 

stretching  frequency  should  be  at  1971  cm  If  the  66  cm  ^  error  for  the  C-0 

mode  is  incorported  into  the  predicted  Rh-H  stretching  frequency,  the  ab 

initio  calculations  suggest  that  an  infrared  band  should  occur  near  1709  am  1. 

This  band  may  be  under  the  wing  of  the  intense  1615  cm  1  band  (spectrum  a) 

which  is  due  primarily  to  the  hydrated  support.  It  is  very  interesting  that 

the  calculations  predict  an  isotopic  shift  of  only  3  cm  1  for  the  C-0 

stretching  mode  upon  deuterium  exchange;  the  observed  value  was  5  cm  ^ 

Further,  the  calculations  predict  that  a  band  should  exist  near  1233  cm  1 

(corrected  from  1169  cm  '*')  corresponding  to  the  Rh-D  stretching  mode  of  the 

carbonyl  deuteride.  There  was  a  prominent  feature  observed  at  1193  cm  1  in 

spectrum  b  of  Figure  2,  but  unfortunately  O-D  groups  on  the  hydrated  support 

would  also  be  expected  to  produce  a  band  in  this  region.  In  fact  when  D20  was 

added  to  the  cell,  a  large  band  peaked  near  1205  an-1.  In  summary,  it  is 

r 

likely  that  we  have  observed  the  Rh-H  band,  but  we  can  not  at  this  point  be 
certain  whether  it  lies  under  the  C-0  band  at  2037  cm-1  or  under  the  hydrated 
support  band  at  1615  cm  Likewise,  the  Rh-D  band  must  either  occur  at  1450 
cm  1  or  under  the  1193/1205  cm  1  O-D  support  band.  Of  course,  the  solution  of 
this  dilemma,  if  and  when  it  can  be  found,  may  determine  the  worth  of  the  ab 
initio  procedure  which  was  employed  for  use  in  future  calculations  related  to 
surface  chemistry. 


Several  other  interesting  points  evolved  from  the  3b  initio  calculations 


performed  in  this  study.  First,  the  ground  state  of  the  rhodium  carbonyl 

hydride  species  (bent  form)  was  predicted  to  be  28.5  kcal/mol  more  stable  than 

the  rhodium  formyl  species.  There  was  a  zero  energy  barrier  predicted  for  the 

conversion  of  the  formyl  species  to  the  carbonyl  hydride  species.  These  facts 

explain  why  the  rhodium  formyl  hydride  species  could  not  be  detected  by  Yates 
14 

and  Cavanagh.  Second,  the  rhodium  carbonyl  hydride  gas-phase  molecule  was 
predicted  to  have  a  bent  ground  state  with  a  HRhC  bond  angle  of  86.6°.  This 
is  gratifying  because  the  carbonyl  hydride  species  on  Rh/Ti02  would  be 
expected  to  be  bent,  and  a  HRhC  bond  angle  near  90°  would  seem  optimum.  Thus, 
we  are  reasonably  confident  that  the  gas-phase  ab  initio  calculations  can  be 
at  least  qualitatively  extrapolated  to  the  supported  catalyst  situation. 

Third,  the  sb  initio  calculations  predict  a  short,  stronger  Rh-C  bond  exists 
for  the  carbonyl  hydride  species  than  for  RhCO.  Thus  there  is  greater 
electron  donation  from  CO  toward  Rh  through  the  Rh-C  O' bond  for  the  former. 
Also,  the  carbonyl  hydride  species  provides  better  overlap  between  Rh  and  the 
if*  orbital  on  CO  because  of  the  short  Rh-C  bond.  These  two  modes  for  electron 
transfer  satisfactorily  explain  why  the  CO  bond  is  considerably  weakened  when 
the  hydride  ligand  is  present.  This  is,  of  course,  the  most  likely  reason  why 
the  presence  of  in  a  reactor  greatly  enhances  the  dissociation  of  a  CO  bond 
for  a  rhodium  carbonyl  species.  It  explains  why  hydrogen  is  necessary  to 
cause  the  dissociation  of  C02  on  Rh, ^  and  it  probably  explains  the  formation 
of  active  carbon  on  Rh  in  the  catalytic  methanation  reaction.  In  some 
experimental  situations  the  dissociation  of  the  CO  bond  in  the  rhodium 
carbonyl  hydride  species  could  be  the  rate  determining  step  for  formation  of 
methane  and  other  hydrocarbons.  It  should  be  noted  that  we  can  not  completely 
rule  out  the  possibility  of  a  transient  formyl  species  being  involved  in  the 


enhanced  CO  dissociation  in  the  presence  of  hydrogen  given  that  the  energy 

barrier  predicted  for  passing  from  the  carbonyl  hydride  species  to  the  formyl 

28 

species  was  only  28.5  kcal/mol. 

Carter  and  Goddard  have  recently  performed  GVB  ab  initio  computations 

29 

concerning  the  insertion  of  into  an  adjacent  Ru-H  bond.  They  found  that 

the  results  were  affected  markedly  depending  upon  whether  electron  correlation 

was  or  was  not  included  in  the  calculations.  At  the  Hartree  Fock  level 

without  Cl  there  was  no  energy  barrier  predicted  for  the  insertion  reaction, 

while  at  the  GVBCI-MGSCF  level  an  activation  barrier  of  11.5  kcal/mol  was 
.  29 

predicted.  We  find  also  that  the  inclusion  of  electron  correlation 
significantly  affects  the  results  for  the  conversion  of  the  rhodium  formyl 
species  into  the  rhodium  carbonyl  hydride  species.  In  fact,  from  Table  II  it 
is  evident  that  the  inclusion  of  Cl  caused  the  energy  of  the  rhodium  formyl 
species  to  increase  by  10.6  kcal/mol.  Also,  an  energy  barrier  was  present  at 
the  HF  level,  but  not  at  the  MP3  level,  for  the  interconversion  between  the 
carbonyl  hydride  and  formyl  species  (see  TS  energies  in  Table  II) .  Thus 
calculations  at  the  HF  level  may  give  misleading  results  for  insertion 
reactions  into  metal  hydrogen  bonds,  and  we  recommend  the  inclusion  of 
electron  correlation  for  such  reactions. 


CrVn 
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The  a£  initio  calculations  in  this  work  have  predicted  that  a  bent 
rhodium  carbonyl  hydride  species  is  much  more  stable  than  a  rhodium  formyl 
species  and  that  there  is  no  energy  barrier  which  would  prevent  the  formyl 
species  from  collapsing  to  the  carbonyl  hydride  species  which  rationalizes  the 
detection  of  the  former,  but  not  the  latter,  by  infrared  for  supported  Rh 
catalysts.  The  calculations  also  predict  that  the  Rh-H  stretching  frequency 
for  the  rhodium  carbonyl  hydride  species  should  occur  some  300  -  400  an  1 
lower  than  the  corresponding  frequencies  for  rhodium  organometallic  complexes . 
Such  a  band  may  exist  for  Rh/Ti02,  but  if  it  does,  it  unfortunately  overlaps 
the  hydrated  support  band,  so  we  can  not  definitely  conclude  that  the  Rh-H 
stretch  does  not  lie  under  the  band  for  the  C-0  stretch  for  the  rhodium 
carbonyl  hydride  species.  The  calculations  do  rationalize  the  ability  of  the 
hydride  ligand  to  enhance  C-0  bond  dissociation  on  supported  rhodium. 
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Tabl**  ii.  Relative  Energies  (kcal/mol)  of  Species  on  the  HRhCO  Potential 
Energy  Surface 


Comparison 

BAS1 

BAS2 

best 

H? 

W2 

IS' 

MP2 

MP3 

^h  -  'W 

47.7 

46.5 

46.1 

46.  lb 

1RhH  -  3RhH 

41.9 

38.1 

36.1 

34.8 

19.7 

21.0 

21.5 

^hCO  -  4RhCO 

21 A 

28.1 

26.6 

29.0 

30.1 

28.4 

28. 4C 

1 

HRhCO 

1.4 

17.8 

14.9 

4.0 

26.2 

17.6 

18.1 

2 

thCO 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

Rh&i 

13.5 

36.4 

26.2 

15.7 

41.6 

26.3 

28.5 

4 

TS 

22.2 

19.4 

18.1 

24.1 

18.9 

20.1 

19.2 

h  -  h 

-55.3 

-62.1 

-65.7 

-49.7 

-49.4 

-54.7 

-54. 7d 

12  -  32 

-56.7 

-79.0 

-71.3 

-55.2 

-91.8 

-73.8 

-73. 8d 

■4  -  33 

29.5 

26.4 

23.6 

22.8 

6.8 

9.9 

9. 9d 

14  -  34 

0.1 

-14.2 

-11.5 

-2.3 

-28.9 

-16.3 

-16.3d 

a)  A  comparison  is  made  of  the  energy  of  the  first  state  relative  to  the 
second  state.  Geometries  on  the  singlet  HRhCO  potential  energy  surface 
(1-4)  are  relative  to  2.  Singlet  -  triplet  splittings  of  1  -  4 

are  computed  at  the  optimized  singlet  geometries.  A  positive  value 
indicates  that  the  higher  multiplicity  is  more  stable. 

b)  A  contaminated  spin  state  is  obtained  for  the  doublet  state  with  BAS2. 
Best  value  is  MP3/BAS1. 

c)  Zero  point  correction  is  not  included.  A  contaminated  state  is  obtained 
for  the  doublet  state  with  BAS1  and  EAS2. 


d)  Triplet  energy  is  computed  at  the  singlet  geometry.  Zero  point  correction  is 
not  included. 


Table  III.  Computed  Bond  Energies  (kcal/mol)  at  Several  Levels  of  Theory 


BAS1 

BAS2 

a 

Bond 

EF 

MP2 

MP3 

HF 

MP2 

mp3 

besta 

3Rh — H 

49.7 

57.9 

60.5 

61.8 

73.5 

74.7 

73.1 

1 

8 

-7.3 

0.1 

2.1 

2.3 

9.5 

11.0 

11. 0b 

XPh — §h 

7.6 

33.5 

37.2 

25.2 

68.2 

63.1 

61.3 

1 

1 

RhCO 

42.8 

80.1 

77.4 

53.9 

112.3 

96.0 

96. 0b 

Rh — CO 

-14.2 

22.2 

19.0 

-5.6 

48.2 

32.3 

28.6 

a)  At  the  MP3/BAS2  +  ZPC/BAS1  level  of  theory. 

b)  Best  Level  is  MP3/BAS2.  Zero  point  correction  is  not  included. 


IS 


a  -1 

Table  IV.  Vibrational  Frequencies  (an  ) , 

Deuterium  Isotope  Shifts 

and  Intensities 

(km/mol)  at  the  BAS1  level 

Molecule 

Type 

Sym 

Freq (Shift) 

Intensity 

Exp 

CO 

C-0 

cr 

2034 

146 

2143b 

HCO 

H-C-0 

a' 

1098(862) 

42 

1090 (852) ° 

C-0 

a' 

1680(1659) 

81 

1861(1800) 

C-H 

a' 

2793  (2026) 

28 

2488(1937)° 

1RhH 

Rh-H 

<r 

1333(948) 

695 

3RhH 

Rh-H 

cr 

1008  (717) 

377 

HRhCO 

Rh-C-0 

ft' 

234  (209) 

160 

Rh-C 

or 

311(310) 

14 

H-Rh-C 

if 

697 (561) 

75 

Rh-H 

cr 

1433(1019) 

2329 

C-0 

cr 

2010(2009) 

644 

KhCO 

Rh-C-0 

a' 

282  (261) 

17 

out-of-plane 

a” 

327 (326) 

2 

Rh-C 

a' 

343(338) 

36 

H-Rh-C 

a’ 

753  (589) 

100 

Rh-H 

a' 

1643(1169) 

1184 

C-0 

a' 

1971(1968) 

692 

2037  (2032) C 

Rh?H 

Rh-C-0 

a' 

167 (165) 

21 

Rh-C 

a' 

352  (324) 

2 

out-of-plane 

a" 

674 (559) 

17 

Rh-C-H 

a ' 

1319(1015) 

52 

C-0 

a' 

1559(1555) 

268 

C-H 

a' 

2604(1905) 

302 

TS 

H-Rh-C 

a' 

577 ie 

488 

Rh-C-0 

a' 

386 

17 

out-of-plane 

a” 

424 

7 

C-H 

a' 

811 

232 

Rh-H 

a' 

1331 

419 

C-0 

a” 

1805 

534 

-1. 


a)  All  frequencies  have  been  scaled  by  0.89  (see  ref.  19);  shift  refers  to 
the  frequency  for  the  analogous  deuterated  species. 

b)  Ref.  20. 

c)  Ref.  21. 

d)  This  work . 

e)  Imaginary  transition  vector  (negative  eigenvalue  of  the  force  constant  matrix) . 


FIGURE  CAPTIONS 

o 

Figure  1.  Geometric  parameters  {bond  lengths  in  A,  bond  angles  in  degrees) 
for  the  species  for  which  computations  are  reported  in  this  work. 

-2 

Figure  2.  Infrared  spectra  of  a  2.2%  Rh/Ti02  film  (4.3  mg  cm  )  following 

—6 

reduction  in  at  480  K,  evacuation  for  1  h  at  490  K,  and  evacuation  to  10 
Torr  at  298  K  overnight  then  (a)  introduction  of  77.4  torr  of  a  1:4  mixture  of 
CQ2/H2,  heating  at  440  K  for  1  h,  and  evacuation  to  10  torr  at  298  K;  (b) 
after  contact  with  100  torr  D2  at  298  K  overnight. 
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